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We show, that photon non-linearities in electromagnetically induced transparency can be at least 
one order of magnitude larger than predicted in all previous approaches. As an application we 
demonstrate that, in this regime they give rise to very strong photon - photon interactions which 
are strong enough to make an experimental realization of a photonic Mott insulator state feasible 
in arrays of coupled ultra high-Q micro-cavities. 
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Introduction: Various quantum mechanical effects 
have been observed with photons. Examples include 
quadrature squeezing and entangled photon states 
All these effects require strong non-linear interactions 
between photons for their observation. Unfortunately, 
photons tend to interact only weakly and it is thus of 
considerable interest to develop schemes achieving larger 
and larger photon non-linearities, as these would make 
further quantum effects accessible in experiments. 

A scheme for the generation of a large photon non- 
linearity has been proposed by Imamoglu and co-workers 
Q for atoms with a level structure considered in elec- 
tromagnetically induced transparency (EIT) which 
interact with the resonant light mode of a cavity and in 
Q it was clarified when this is strictly a photon-photon 
interaction. Non-linearities obtained in this way can be 
strong enough to employ them as single photon turn- 
stile devices [5[- The photon non-linearity considered in 
Q can be studied in a more general regime as a non- 
linear interaction for dark statepolaritons [f| . This non- 
linearity and other approaches [7| have been employed to 
show that effective many body physics should be observ- 
able in arrays of coupled cavities [6[. 

In this letter, we show, that the photon non-linearity 
proposed in Q exists in a much more general parameter 
regime than previously considered [H]. In that way, we 
can relax a restricting assumption used in 0, H| and Q 
and are thus able to increase achievable photon-photon 
interactions by at least one order of magnitude. As the 
non-linear interactions presented here are significantly 
stronger than in previous schemes they open up possibil- 
ities for the observation of phenomena which were pre- 
viously not accessible in experiments. If generated in an 
array of coupled ultra high-Q micro-cavities [1], they can 
become large enough to make the experimental observa- 
tion of a Mott insulator state for photons rather than 
polaritons feasible. 

The parameter regime we consider in this paper is cru- 
cially different to the dark state polariton regime em- 
ployed in Q and the regime addressed in Q. Indeed, the 
non-linearity we derive here is an interaction between 
photons and not polaritons (joint atomic - photonic ex- 



citations) as considered in 0, 0, @, 0] • 

In the photonic Mott state we consider, exactly one 
photon exists in each cavity, provided the whole structure 
contains on average one photon per cavity || . Moreover, 
the photons are localized in the cavity they are in and 
are not able to hop between different cavities. In such a 
situation photons behave as strongly correlated particles 
that are each "frozen" to their lattice site, a system that 
would correspond to a crystal formed by light. The non- 
linearity derived here, when applied to coupled cavities, 
could thus provide the first realistic possibility to observe 
light in this exotic quantum state. 

As this photonic Mott insulator is characterized by the 
number of photons in a single cavity being, (ajai) ~ 1 
for cavity /, and its fluctuations being zero, ((aja;) 2 ) — 
(ajai) 2 ~ 0, it is fundamentally different to polaritonic 
Mott insulators as predicted in @ and 0, where the 
number of polaritons at one site is exactly unity but the 
photon number fluctuates ({ajai) « 5 and ((ajai) 2 ) — 
(ajai) 2 i for no detuning). 

This letter provides two main results: We general- 
ize the photon-photon interaction 0] showing that much 
stronger interactions can be obtained and, as an example 
for the usefulness of this improvement, we demonstrate 
that with this new tool, the experimental observation of 
a Mott insulator made of photons/light is within reach 
with present day high-Q micro-cavities. 

Photon-photon interaction: We begin by considering 
one cavity that interacts with 4-level atoms which are 
driven by an external laser in the same manner as in 
EIT, see fig. Q] The transitions between levels 2 and 3 are 
coupled to the laser and the transitions between levels 2- 
4 and 1-3 couple via dipolc moments to the cavity mode. 
The existence of a strongnon-linearity in this atom cavity 
system was shown in @, 13] and a similar nonlinearity has 
recently been observed experimentally Q . 

In a rotating frame with respect to Hq = oj(a^a + 

■i) + YljLi ( UJ(J 22 + Wf7 33 + 2wcr^ 4 ), where we set h = 1 
throughout the paper, the Hamiltonian of the atoms in 
the cavity reads, Hj = J2f=i i £<7 22 + ^°33 + (A + s)a{ 4 ) + 
Ej=i(^°23 +513 013 flt +524 024 flt + h.c.) , where ct^ = 
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FIG. 1: The level structure and the possible transitions of 
one atom, uj is the frequency of the cavity mode, fi is the 
Rabi frequency of the driving by the laser, $13 and <?24 are 
the parameters of the respective dipole couplings and S, A 
and e are detunings. 



\kj)(lj\ projects level / of atom j to level k of the same 
atom, lo is the frequency of the cavity mode, Q is the 
Rabi frequency of the driving laser and 513 and 524 are 
the couplings of the cavity mode to the respective atomic 
transitions which are all assumed to be real. S, A and 
e are detunings of atomic transitions with respect to the 
cavity and laser fields, see fig. [1] All atoms interact 
in the same way with the cavity mode and hence the 
only relevant states are symmetric Dicke type dressed 
states 

When the system is operated in the strong coupling 
regime, atoms and photons do no longer behave as sep- 
arate entities and polaritons, the excitations of dressed 
states, become a useful and intuitive picture for its de- 
scription. For the derivation, we thus make use of 
the polariton picture considered in [|| but will later 
consider the photonic limit. Hence while the tech- 
nique is similar the physical conditions are different to 
those considered in Q. We define the polariton cre- 
ation operators Q, p\ = B^ 1 \gS\ 2 — f2au and p± = 



A(A±S) 



where g = vNgiz, 
B = v^TTF, A = + S 2 , S\ 2 = ^Ef=i<4 

and S\ 3 = ^2j=i °3i • The operator p\ excites the 
dark state in which level 3 is not occupied because both 
ways to excite it interfere destructively Q . Note in par- 
ticular that in the limit Q 3> g, the dark state polaritons 
Pq become photons. 

Neglecting the coupling to level 4, g 2 4 = 0, and two 
photon detuning, e = 0, the Hamiltonian Hj in the rel- 
evant subspace spanned by Dicke type dressed states is 
given by [Hj] 0tE=Q = fi plpo + H+P+P+ + M- P-P- , 
with the polariton frequencies /io = and fi± = (6±A)/2 
Hence Pq,p\. and p\_ decouple for £ = g 2 4 = 0. 

In a regime where A, 1724 <C £1, only the dark state 
polaritons p^ couple to the atomic level 4 as for the 
other species their separation in energy from level 4 is 
much larger than their coupling to it. Hence the cou- 



24 t h.c. can be approximated to read 

— 524 g ^ B~ 2 S\ 4 Pq + h.c. @ . In the limit g <C fl we can 
furthermore expand the coupling in powers of g to obtain 
— fl24fl^ _1 S\ A a 2 + h.c. Provided that fl24 fl <C |Afi|, 
this coupling can be treated perturbatively, yielding a 
2nd order energy shift of n p h (n p h — 1) g^g 2 A -1 , 
where n p h is the number of photons in the cavity. Note 
in particular that the assumption 024 <C A, which has 
been used in earlier approaches [2|, [6|, is not made here. 
This assumption is not necessary since the atomic level 
2 has very small occupation in the limit !1 > 5 and the 
perturbative treatment is justified even if fl24 > A. This 
gives rise to photon non-linearities that can be orders of 
magnitude larger than in previous approaches 0, @1 ■ The 
photon non-linearity can be expressed as the effective 
Hamiltonian — fl| 4 g 2 A" 1 Vl~ 2 a' a' a a . Note in partic- 
ular that this photon-photon interaction is repulsive for 
A < and attractive for A > 0, giving easy experimental 
access to both regimes. An analogous derivation shows 
that the two photon detuning e (see fig. [T]) leads to a shift 
of the photon frequency given by the effective Hamilto- 
nian £ g 2 il~ 2 cr a, which acts as a chemical potential. 

The more the dark state polaritons become truly pho- 
tons (g -C fi), the larger is the gain in non-linearity com- 
pared to previously considered scenarios @, @1- If we 
assume 51 = lOOfl and |A| = 524/10, the requirements 
g < fl and fl 2 4 g / | A fi| = 1/10 <C 1 hold. If one followed 
previous derivations one would need |A| fl 2 4 and 

therefore choose |A| ~ \Qg24 or larger. Since the non- 
linearity is proportional to 1/A, it is in our approach a 
factor 100 larger than in previous ones. 

An effective decay rate T for the photons can be ob- 
tained following a similar derivation as in Q and only 
keeping leading terms in g/fl: T = Tc + 6(n p h — 
2) g\ A g 2 1ST 2 ?T 2 T 4 , where 0(n ph — 2) is the Heavi- 
side step function, Tq the bare cavity decay rate and 
T4 the rate for spontaneous emission from level 4, which 
only occurs if 2 or more photons are present in a cavity 
and is further suppressed since fl24fl <C |AQ|. Cavity de- 
cay is thus the main source of photon loss Hence, 
although the non-linearity achieved here is much larger 
than in @, Q , it does not come at the expense of stronger 
decay rates. 

Photonic Mott insulator: We now consider a periodic 
array of cavities. The resonant modes of neighboring 
cavities have a nonvanishing overlap which gives rise to 
photon hopping between them [l]| . Including the on-site 
potential, the complete Hamiltonian for e = takes on 
the form of a Bose-Hubbard (BH) Hamiltonian. 
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Here Y1(r r>) * s t ne sum of all pairs of cavities which are 
nearest neighbors of each other, J is given by an over- 
lap integral over neighboring cavity modes and can be 
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obtained numerically for specific models, creates a 

photon in the cavity at site R and U = <?f 4 g 2 A -1 fi~ 2 is 
the on-site potential. For successfully observing a Mott 
insulator state, the repulsion term U needs to be much 
larger than the damping T. This can for example be the 
case for toroidal micro-cavities @ , for which the achiev- 
able parameters are predicted to be c/24 ~ 2.5 x 10 9 s _1 , 
T 4 - 1.6 x 10 7 s- 1 and T c ~ 0.4 x 10 5 s _1 . Taking 
g 24 g/\An\ = 0.1 and g/ft = 0.1, we obtain U/T ~ 625. 

To provide evidence for the accuracy of our approach 
and to show that the repulsion U can indeed be strong 
enough to observe a Mott insulator state for photons, we 
present a numerical simulation of the dynamics of 3 pho- 
tons in 3 coupled toroidal micro-cavities that can tunnel 
between neighboring cavities according to J(a\ a%, + h.c.) 
(a\ creates a photon in cavity 1 etc.) and interact 
with atoms in each cavity. We take parameters, which 
are predicted to be achievable in and assume S = 
1.0 x 10 u s-\ A = -1.25 x 10 9 s-\ ,913 = 2.5 x lOV 1 , 



N = 1000, e = 0, J = 1.2 X lO^- 1 and fi = 20VNgi 3 = 
1.58 x 10 12 s-\ which yields U = 1.24 x lOV 1 . The ini- 
tial state is 1 photon in each cavity, i.e. I^o) = |1, 1, !)■ 
Cavity decay and spontaneous emission are taken into 
account in the standard way [l|. We compare this dy- 
namics to that of 3 particles in the 3-site version of 
(fl]) with additional decay terms at rate T focussing on 
the expectation value and the fluctuations of the num- 
ber of photons in one cavity, i.e. ni = (ajai) and 

Ft = ^((a^) 2 )-^) 2 (I = 1,2,3). 

Fig. [2] shows the dynamics of 3 photons in 3 cavi- 
ties with periodic boundary conditions. Figures [2] a and 
b show the number of photons, m, and number fluc- 
tuation, Fi, in one cavity according to the full model. 
Figure [2] c shows the differences between full and effec- 
tive BH-model, 5m = ["-i] cavitics — ["i] BH (solid line) and 
SF t = [Fi] cavitics - [Fi] BH (dashed line). All these three 
plots show results for one single trajectory of a quantum 
jump simulation (Til ] . The abatement of the norm showed 
that the damping rate T accurately describes the model's 
losses (deviations less than 5%). Figure[2]d finally shows 
ri\ and F\ as given by a solution of a master equation 
for the 3-site BH-model including damping according to 
r, which corresponds to the values to be observed in an 
experiment. There is a probability of ~ 8% to loose the 
photon in the simulated time range. 

Mott insulator to superfluid transition: Due to the 
high Q-factors of present day micro-cavities our approach 
even allows to observe quantum phase transitions for in- 
teracting photons. We consider here the transition from 
the Mott insulator state to a superfluid regime. The on- 
site potential U can be decreased by increasing the inten- 
sity of the driving laser and hence fi. Fig. [3] shows that 
the local photon number fluctuations are very small in 
the Mott phase and increase as U becomes smaller than 
J and the system is driven into the superfluid phase. 






-0.01 



-0.03 
-0.04 





11 1; ji w 


^ViT 1 1 1 
" 1 


I \ 1 1 1 1 
1 11 1 , , 



0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

f in s t in 10~ 6 s 

FIG. 2: Mott insulator state for three polaritons in three 
cavities: a: photon number n\ , b: photon number fluctuation 
Fi, c: differences between full and effective BH-model, 5m 
(solid line) and 8F1 (dashed line) and d: rzi and Fi as given 
by a solution of a master equation for the three site BH-model 
p| including damping V. 
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FIG. 3: Mott insulator to superfluid transition for three po- 
laritons in three cavities: a: On-site potential U (solid line) 
and hopping J (dashed line), b: number of photons in cavity 
1, 711 for the full model (solid line) and the effective BH- model 
(dashed line), c: fluctuations of the photon number in cavity 
1, Fi for the full model (solid line) and the effective BH-model 
(dashed line) and d: ni and Fi as given by a solution of a 
master equation for the three site BH-model ([T]) including 
damping F. 



The parameters are the same as in figure [5] except for 
A = -2.5 x lOV 1 and J = 2.5 x lOV 1 . fi is initially 
fi, = 10v7Vgi3 and is increased to fi/ = 100%/iV 513 to 
drive the system through the transition. 

Let us stress again, that figures [2] and [3] show numbers 
and number fluctuations for photons and not polaritons. 
The photon number statistics in one cavity can be mea- 
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sured via a procedure described in [l5| and p] (c.f. [16[ 
for applications). 

Experimental requirements To generate a non- 
linearity that sufficiently exceeds the decay rates such 
that the resulting Mott state has a long lifetime, atom- 
photon couplings should exceed cavity decay rates by one 
or two orders of magnitude. Promising candidates for 
an implementation are therefore toroidal micro-cavities 
which have a very large Q-factor (> 10 s ), can be pro- 
duced in large numbers, positioned with high precision 
and efficiently coupled to optical fibres as well as to Cs- 
atoms via their evanescent field Q. Photonic crystals 
represent an appealing longer term alternative as they 
offer the possibility for the fabrication of large arrays of 
cavities in lattices or networks |lOI| . Other recently devel- 
oped micro-cavities, such as [12j [. have very high coopcr- 
ativity factors but are less suited for the present purpose 
because of their relatively high cavity decay rates. 

In contrast to Q and there is no restriction on N as 
excitations are photonic and the requirement that level 2 
is long lived can be relaxed as it is very weakly occupied. 
Therefore, circuit cavities interacting with a Cooper pair 
box 11| which would have very low cavity decay if op- 
erated in a row (Q ~ 10 6 ) are very suitable for an im- 
plementation, too. Here, a ratio of U/T ~ 15 could be 
achieved in current experimental devices. Current mea- 
surement schemes in circuit cavities would however re- 
quire open end ports implying higher cavity decay. 

Variations in 313 or N for different cavities can be com- 
pensated by individually adjusting the laser driving f2 for 
each cavity. Relative fluctuations in ^24 of as much as 
50% are permitable as U cx 324 and the U of each cavity 
would still be much larger than T but much lower than 
\/j,± — /i |, (\fx± — /io| ~ 10 12 s _1 for the parameters of fig- 
ure [2]). For the same reasons relative fluctuations of 50% 
in A are tolerable. A possible detuning 5c between neigh- 
boring cavities reduces the photon hopping by a factor 
J /5c, thus even favoring the Mott state, and is harm- 
less if the energy separation between polariton species 
persists, i.e. 5c <§C \n± — Mo I (8c < 10 10 s _1 for micro- 
toroids). However only for \U\ > \8c\, the non-linearity 
U is the only origin of the immobility of the photons. De- 
tunings from the two photon resonances with the driving 
laser generate effective chemical potentials, see above. 
Variations in J are tolerable as long as J <C |/i± — Mo I- 
We note that signatures of the photonic Mott state, such 
as Fi f» in contrast to the superfluid regime, can already 
be observed for as little as two coupled cavities. 

The demonstration of the strong coupling regime in 
micro-toroidal cavities 0] was done with atoms falling 
through the interaction region, where the interaction 
times were about 2/xs. For an experiment these cavities 
could thus be loaded with photons forming a photonic 
superfluid that are stored long enough due to the high 
Q. The atoms are then dropped along the cavities. As 
they approach the evanescent field, their interaction with 



it increases and the system is driven into the Mott state 
where it could remain for the duration of the interaction, 
i.e. 2/is which is enough for our proposal. 

Conclusions: We have shown that photon non- 
linearities in EIT systems can be much larger than pre- 
viously thought. Our approach does not use additional 
technology such as lasers etc but demonstrates that cer- 
tain constraints which were previously thought to be nec- 
essary may be relaxed. With the non-linearities we pre- 
dict, Mott insulator states and the quantum phase tran- 
sition between superfluid and Mott insulator become ob- 
servable for photons trapped in high-Q micro-cavities. 
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